function u_cara = u_cara(gamma,x)
u_cara = -1/gamma*exp(-gamma*x);
end
